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^ ■ Abstract 

g . We investigate diffusion-limited aggregation (DLA) in a wedge geometry. Arneodo and collab- 

\ orators have suggested tfiat the ensemble average of DLA cluster density should be close to the 

O ' noise- free selected Saffman- Taylor finger. We show that a different, but related, ensemble average, 

o : 

that of the conformal maps associated with random clusters, yields a non-trivial shape which is 
^ , also not far from the Saffman- Taylor finger. However, we have previously demonstrated that the 

' same average of DLA in a channel geometry is not the Saffman- Taylor finger. This casts doubt on 

in ■ 

the idea that the average of noisy diffusion-limited growth is governed by a simple transcription of 

O , 

'lO ■ noise-free results. 
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The diffusion-limited aggregation model uses aggregating random walkers to 
form a random fractal pattern. We can also use random walkers to generate a 
conformal map from the exterior of the cluster to the exterior of the unit circle. 
The map contains all the information about the growth probabilities for various 
points on the cluster. We define the pattern generated by the ensemble average 
of the conformal maps as a kind of mean-field pattern for DLA. In previous 
work we showed that there seemed to be an intriguing relationship between 
DLA in a channel and the noise-free Saffman- Taylor finger. We show here that 
this relationship is much more ambiguous when we apply it to DLA in a wedge 
geometry. This result casts doubt on the 'averaging conjecture' which holds that 
the average of noisy growth 'remembers' noise-free results. 



I. INTRODUCTION 



The last century has abounded with examples of unexpected richness in the problem of 
brownian motion as formulated by Einstein in his seminal paper in 1905 [ll. In the same 
year, Pearson [2| pictured the process as a random walk. In the past century random walks 
and brownian motion have become central themes of statistical physics. The problem is 
particularly astonishing in that it constantly generates new ways to think about nature, 
and new descriptions of physical processes. One example was the discovery by Witten and 
Sander ^, Q| that aggregating random walkers give rise to random fractal patterns. This 
process, diffusion-limited aggregation (DLA) is simple to formulate: a seed particle is put 
at a point, and then a random walking particle is launched and allowed to proceed until 
it touches the seed; then it stops. Then another walker diffuses until it touches one of the 
first two, and so on. We study the cluster generated. This is the simplest paradigm of 
noise-dominated growth limited by diffusion, a common natural process. 

This process has been studied intensively over the last twenty-three years (for a review, see 
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but there is still not complete theoretical understanding, though some recent progress 

n 

has been made [6] in describing local correlations. In this paper we investigate a property 
of the ensemble of DLA clusters namely the generation of average shapes. Our motivation 
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comes from the remarkable suggestion of Arneodo and collaborators 0,0] that the (somewhat 
arbitrarily defined) ensemble average shape of noisy DLA clusters, would be the pattern 
generated by noise-free diffusion-limited growth. We will refer to this as the averaging 
conjecture. The relevant case of noise-free growth is the Saffman- Taylor viscous finger lollioll. 
i.e, the shape of the surface of an inviscid fiuid invading a viscous one, as water into oil. 
The work of Arneodo, et al. showed that Saffman- Taylor fingers in a channel and a wedge 
were close, but not exactly the same, as DLA averages. However the work had a number of 
arbitrary parameters, and, in the wedge, there were serious ambiguities. 

We investigated the averaging conjecture in our work ll| on DLA in a channel with 
refiecting boundaries at the walls. There, we formulated a new definition of the average 



shape by averaging the conformal map 



13l 1^ that generates the cluster - this amounts 



to weighting points on the surface according to their growth probability. In that work we 
found that the average shape of DLA (using our definition) was not a Saffman- Taylor viscous 

"hL., a ..... . a... B Q .... ..... .... ... ..o. 

fingers are actually closely related growth processes; DLA does not have surface tension like 
viscous fingering, but rather a fixed particle size that defines the tip radius. We followed 
up the suggestion 0, Q that fiuid fiow with surface tension is closely related to a variant 
of DLA called the dielectric breakdown model (DBM) 17| with parameter r/ ^ 1.2 (to be 
defined below). In fact, the suitably averaged DBM clusters then turned out to fit the 
Saffman- Taylor shape quite closely. 

Here we look at the related problem of DLA in a wedge. Previously jilSi] we have shown 
that tip-splitting in this geometry gives us access to local correlations. In the present context 
of looking at average shapes, the wedge is interesting in several respects. This is a richer 



problem than that of the channel in that the wedge an 



gle is a free parameter, and the shape 



of the fluid invasion is more complex l21l . l22| |. We pose the question of whether 

the DBM shape in a wedge is also a good approximation to the Saffman- Taylor flnger in a 
wedge. 

This study has a special signiflcance in this focus issue on brownian motion. Our work 
here is primarily numerical. The simulations are all done by random walker sampling: we 
generate DLA clusters, DBM clusters, and conformal maps by this single method which, as 
it happens, is by far the most efficient method available. This century-old technique has not 
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FIG. 1: A DLA cluster in a 90° wedge with reflecting boundaries. The structure of the hierarchical 
maps is also shown. 

lost its freshness and power. 

II. DLA AND DBM IN A WEDGE 

There are now available very sophisticated schemes for generating DLA clusters. The one 
we use is based on the method of hierarchical maps js^. In this method space is divided into 
regions of various sizes which help keep track of the nearest points on the cluster. Then the 
random walker can make large jumps in empty regions, vastly speeding up the computation. 
With this method the time to create an particle cluster is proportional to A^^ with p ^ 1.1. 

We need to make DLA clusters in a wedge with reflecting boundary conditions. We 
do this by means of a trick. It there is a wall at some position, every time we deposit a 
walker, we deposit an image walker reflected in the wall. Suppose we are interested in 90° 
wedges. Then we have two perpendicular walls, and four walkers are deposited at once. 
Using this method we can use a radial DLA code to produce wedges of opening angles 
180° /n,n = 1,2, .... In this paper we will concentrate on 90° and 60° wedges. An example 
of a 90° wedge is shown in Figured 

In the following we will see that we need to consider DBM clusters. These are defined as 
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(2) 



follows: we imagine that the cluster is a grounded conductor with unit charge. Outside the 
cluster define a potential, (p so that: 

= 0, 0, = 0. (1) 

This defines an 'electric field' on the surface, dip/dns- Then the dielectric breakdown model 
takes the growth velocity on the surface of the cluster to be: 

dn 

We interpret this equation probabilistically: f„ is taken as the density of the growth proba- 
bility on the surface, /i. In practice, add a particle at a point on the cluster with probability 
proportional to f„. It is known that at large scales DBM clusters with r] = 1 have the same 
scaling as DLA. 

The original method 0| to grow DBM clusters was to solve Eq. by relaxation. This 
algorithm is very slow, and is not practical for generating large clusters. Recently ^24j we 
have introduced a method of growing DBM clusters by random walker sampling. The key to 
this method is to define the age, ai of a growth site. This is the number of random walkers 
that have landed anywhere on the cluster since the last particle grew at the site. We can 
also define a^, the number of walker that have landed since the k^^ most recent particle grew 
there. 

Since the frequency of landing of random walkers is proportional to /x, it is clear that 
1/ai, 2/a2, .. at a site are estimates of /i at that site. We have shown [24] that this estimate 
is adequate to allow us to grow DBM clusters. In our work here we use to estimate the 
probability. 

The method of growth is as follows: if a particle lands at a site with low probability, we 
arrange to have it add little to the cluster, and at high probability sites, add a good deal. 
This is accomplished by adding a mass, 5m = A^^~^ oc a]^^ . (The power is 77 — 1 because we 
already have a probability yU for the walker to land at the site.) In practice, when a particle 
is added at a site, it is moved onto the existing particle so that a portion proportional to 
a]~^ contributes to new growth. We also change the prefactor, A, as we go along to make 



an efficient code. For details, see [2^. Examples of this is shown in Fig. 0. 

The computations we discuss below are averages over an ensemble of off-lattice DLA and 
DBM clusters grown in this way. Our DLA clusters had 1,000,000 particles in the wedge 



FIG. 2: A portion of DBM clusters for two different values of rj. 
particles and the enhanced growth at the tips is easy to see. 



For rj = 2 the overlapping of 



and we averaged over an ensemble of 400 realizations. For the DBM clusters we had 500, 000 
particles in the 60° wedge, 1,000,000 particles in the 90° wedge and also 400 realizations. 
Our motivation for going to these large sizes is the fact that DLA suffers slow crossovers 



25|. If we use small clusters we are not seeing the asymptotic behavior. 



III. CONFORMAL MAPS 

In recent years a number of groups have looked at DLA in an entirely new way. 

The cluster shape is considered to be a grounded conductor, as above, and the complex 
potential, 0, is sought, as in Eq. ((l}. The technique introduced was to define a conformal 
map from the exterior of the cluster to the exterior of the unit circle. The Laplace equation, 
Eq. ((H) can be solved easily outside the circle, and the solution mapped back to the cluster. 
Of course, the solution for d(f)/dn is uniform on the circle. Thus the image of two parts of 
the cluster perimeter with the same growth probability will map to parts of the unit circle 
with the same length. Put another way, the inverse images of uniformly spaced points on 
the unit circle are distributed on the cluster with densit y u . 

In order to construct the map Hastings and Levitov invented an iterative technique 
which grows a cluster and calculates the map at the same time. This is a practical method, 
but slow. In Ref. we constructed an alternative method which is much faster. We grow a 
cluster by the conventional fast scheme using random walkers. Then we freeze the cluster at 
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FIG. 3: The averaged profile of DLA and DBM clusters with rj = 1.2 in a 90° wedge. Also shown 
is the the analytic solution for a Saffman- Taylor finger with the selected A from [igl. [2ol|. 



the desired size, launch n random walkers as probes, and record where they land. Then the 
values of the map on the unit circle are found as follows. We choose one of landing positions 
as a starting point, say along the x-axis, and number the walkers around the perimeter of 
the cluster. By the observation above, if we are at walker m, we know that we have turned 
an angle 6 = 2Tcm/n + 0{n^^^'^) from the image of the x-axis on the unit circle. In Ref. 
jl^ we used this information, which samples the boundary values of the conformal map, to 
construct the map itself by analytic continuation. 

For our purposes here we use the information in a simpler way. We constructed the map 
to the unit circle using n = 100, 000 for each cluster that we grew. For each member of 
the ensemble there is a point, r{6) whose image is a point on the unit circle at e*^. Our 
definition of the ensemble average shape generated by the DLA or DBM process j^l is the 
ensemble average of r{9), i.e., the centroid of those points. In Figures (jH)) and (jlj), below we 
show the average shapes in 90° and 60'^irc wedges. 

This definition has a number of advantages: it is unambiguous, in contrast to definitions 
based on density profiles (see below). It is an average which is weighted by probability. That 
is, we are sampling where the growing tips of the cluster are located. It might be accessible 
to theoretical investigation since the growth process is the defining property of DLA. 
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IV. THE AVERAGING CONJECTURE AND MEAN-FIELD THEORY 



A. Saff man- Taylor fingers 

The investigation of Hele-Shaw flow in a channel is an old subject and a good deal is 
known about it. In particular, the invasion of a viscous fluid by an inviscid one is one of the 
most famous of all pattern formation problems 0, . The general result is that the inviscid 
fluid forms a finger, the Saff man- Taylor finger, which fills a fraction. A, of the channel. For 
each A in [0, 1] there is a solution to the fiuid-fiow equations without surface tension. For 
small surface tension, the pattern is very close to that with zero surface tension and A = 1/2. 
This is called the selected value. The finger elongates in time, and is of constant shape in a 
moving reference frame. 

The theory of this sort of flow is remarkably well developed. A few basic notions are 
necessary here: the fluid velocity, v, is derivable from a potential in this case, and is governed 
by D'Arcy's law. The upshot of these two facts is that: 

V • V = = V^0; oc 

That is, fluid flow obeys the same equations as DBM (cf., Eq. (0)) with rj = 1 |26|, so that 
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the Saffman- Taylor finger has a good deal in common with DLA. The bubble of inviscid 
fluid plays the role of the aggregate. However, there are two salient differences: DLA does 
not have surface tension but, rather, a finite size cutoff; and DLA is dominated by noise, 
whereas the Saffman- Taylor finger is a stable, noise-free pattern that is observed for slow 
flows in a channel. 

There is another problem related to the channel problem which also admits of an exact 
solution, that of viscous fingering in a wedge. In this case, for short times there is also a 
selected shape in experiments, at least for a finite time Tu and Ben Amar 20, 21, 2^ 



worked out the theory in this case and showed that there is a self-similar shape that is 
selected. It is a non-trivial pattern whose form is given by a differential equation that needs 
to be solved numerically. Associated with the problem is a selected angle, defined as the 
opening angle at the base of the wedge; see Figures Q and (jD). The ratio of the opening 
angle of the finger to that of the wedge is also called A. With no surface tension there is 
a solution for all A, but, once more, there is a selected value which depends on the wedge 
angle. There is a complication in the fiuid-fiow problem. If the inviscid fluid is pumped 
in a constant pressure, as time goes on there will always come a point when the pattern is 
unstable against tip-splitting. 

B. Averaging 

Arneodo and collaborators 0, 0] exploited the resemblance between DLA growth and 
viscous fingering in the following way: they speculated that the average of many DLA 
clusters would, in some sense, remove the noise, and recover the noise-free pattern. Since 
DLA has no surface tension, but rather a fixed particle size (playing, roughly, the role of the 
capillary length) they assumed that the limit of small surface tension was the appropriate 
one. This is what we refer to as the averaging conjecture. 

They tested the conjecture by generating on-lattice DLA clusters in a channel and aver- 
aging the density, point by point |7[. The density average, p(r) is a function that goes to 
zero at the edge of the channel, and has a maximum, Pmax at the center. One of the level 
sets of this function, Fg : {p(rs) = \pmax} traced out a Saffman- Taylor finger with A = 1/2. 

However, closer scrutiny made the picture more complex. Lattice effects are known to 
distort DLA clusters, and are irrelevant to the kind of physics being considered. Therefore, 
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they returned to the problem j8| and generated off-lattice DLA clusters. They found that 
the level set at \pmax filled 56% of the channel rather than 50%. Or, alternatively, the level 
set that was needed to make A = 1/2 was at Q.Qpmax- 

For the same problem in wedges of various opening angle, on-lattice DLA density averages 
once more gave remarkable agreement with selected fingers j.7']. However, off-lattice the 
situation was different: tip-splitting made the front of the fingers too fiat, though there was 
qualitative agreement with the opening angle of the average density. And, it was necessary 
to choose a level set somewhat arbitrarily. In fact, since the overall density decreased as 
the length of the finger increased, it was necessary to define the opening angle by taking a 
fraction of the maximum density at that distance from the apex of the wedge. 



C. Mean- field theory 



revised the 



Inspired by the averaging conjecture, Levine, Tu, and collaborators 
mean-field theory of Witten and Sander to attempt to write a proposed set of equations 
for the mean density of a cluster. They found qualitative agreement with Ref. In Ref. ^ 
the theory was extended, and other work has followed this up more recently 29|. 

The salient features of this theory are that a level set, defined as above, more-or-less fits 
the channel finger, but that the wedge-based fingers are too fiat. 



D. Average conformal maps for DLA in a wedge 

The numerical underpinnings of the averaging conjecture are troublesome in several ways. 
At the most simple level, the number of particles in the clusters studied were very small, of 
order 10^ — 10^. We have already remarked that DLA in that regime is far from asymptotic. 
More significantly, there are far too many fitting parameters in the discussion. The level-set 
is chosen arbitrarily, and, for many of the discussions, A is chosen to fit the pattern. 

We returned to this problem with our new definition of the ensemble average pattern, 
described above. In a channel we were able to show that the average shape does not fit 
any Saffman- Taylor finger. The finger width of the pattern corresponded to A ~ 0.6, just 
as in Ref. Q], despite the different definitions of the average. With the resolution that we 
had available, we were able to show definitively that the finger we generated did not fit the 
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Saffman- Taylor pattern for any A. 

Here we return to the problem for growth in a wedge. Using the techniques described 
above, we have grown DLA clusters and averaged for wedge angles of 90° and 60°. The 
results are shown in Figures Q and (@)) along with the analytic solutions for the selected 
finger shape. Now, quite remarkably, the opening angle of the finger does fit rather well to 
the selected A. Note that there are no adjustable parameters in this fit whatsoever. The 
tip of the finger in the 90° wedge is flatter than it should be. However, for 60° the fit is 
reasonably good. 



E. Average conformal maps for DBM in a wedge 

In our channel work we tried to salvage the averaging conjecture in the following way. 
There is a theoretical suggestion 0, Q| that the correct analogy between fluid flow and flow 
with surface tension was not with DLA. Rather, the nature of the short-range cut-off in 
Hele-Shaw flow is that the tip radius, R obeys the relation Rv^^"^ = const, rather than the 
DLA case of Rv*^ = const. This led to a relationship between models with different cutoffs 
and different rj such that the dominant growth probability was the same. The model with 
fixed R that grows in the same manner as in fluid flow corresponds to ?7 ~ 1.2. We were 
very encouraged in the channel by this analogy because the average of DBM with this rj did 
fit a Saffman- Taylor finger with A = 1/2 rather well. 

In the wedge, as Figures and (j3)) show, the DBM fingers do not fit the analytic shape. 
They tend to be too narrow. This is exactly the opposite of the situation in the channel. 
Once more, there are no adjustable parameters available to us to fit. 



V. SUMMARY AND DISCUSSION 



When we began this study we were confident that we would be able to put the averaging 
conjecture on a firm footing, based on our experience with the channel geometry. Our 
expectations were not at all fulfilled. It is possible to maintain that we have, in fact, 
eliminated the averaging conjecture altogether. Perhaps this strong conclusion is premature, 
but, certainly, the situation is not very clear. To summarize: in a channel, averaged DBM 



using the mapping of Ref. 



16| gives a Saffman- Taylor finger with the correct A but 
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averaged DLA does not. In the wedge, averaged DLA gives the correct finger opening angle 
for both wedges that we looked at, but averaged DBM does not. 

We might be tempted to say that our proposal for averaging based on conformal maps 
should simply be discarded. If we do that, we are reduced to using density averages which 
don't fit the analytic fingers any better (worse, in fact) and are ambiguous to boot. 

In passing, we should comment on the situation with tip-splitting. We do not agree that 
DLA averages in a wedge should tip-split for any opening angle > 0. We base this on our 
work in the wedge geometry [3| • In that paper we used not refiecting boundary conditions, 
but periodic boundaries for the wedge. We looked at angular correlations of the density for 
DLA clusters, and found that for small wedge angles there was a minimum in the correlation 
function half-way between the branch and its image. We interpreted this by saying that for 
small angles there was one major branch. For large angles, > 144° we found a secondary 
maximum in the correlation function, i.e. more than one major branch. These results do 
not directly carry over to the present case, but, qualitatively, we think that tip-splitting in 
the sense lust described is not at all clearly present in this case. There is a numerical result 
of Ref. js| for a 60° wedge which seems to contradict this, but we are confident that our 
statistics were much better. 

It is likely that tip-splitting for DLA is a probabilistic matter. We suspect that in a 
90° wedge some clusters split, but the majority do not. This could account for the small, 
but definite disagreement between the shape of the tip of the averaged clusters in the 90° 
wedge with the analytic solution. However, there is almost certainly a qualitative difference 
between DLA and fiuid fiow with respect to tip-splitting. Also, the extent of tip-splitting 
found in mean- field theory ^, 28| does not agree with our DLA averages, or with the results 
of 0. 

The project of formulating a description of the average over the DLA ensemble still seems 
to us to be quite a worthy one. However, the present results show that the current state of 
the art in this area is far from giving the definitive answer. 
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